LONG-TIME ASYMPTOTICS FOR THE TODA LATTICE IN 
THE SOLITON REGION 
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O 

o 

^vj ' Abstract. We apply the method of nonUnear steepest descent to compute 

' the long-time asymptotics of the Toda lattice for decaying initial data in the 

^~i^ soliton region. In addition, we point out how to reduce the problem in the 

^ , remaining region to the known case without solitons. 
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1. Introduction 

In this paper we want to compute the long time asymptotics for the doubly 
infinite Toda lattice which reads in Flaschka's variables (see e.g. [16], [l7], or [18]) 



_H , .^^. bin,t) = 2{a{n,t)^~a{n~l,t)^), 

^,' d{n,t) =a{n,t){b{n + l,t)-b{n,t)), 



(n, i) G Z X K. Here the dot denotes differentiation with respect to time. We will 
^ , consider solutions (a, b) satisfying 

a^. (1-2) V|n|'(|a(n,t)--| + |6(n,t)|)<oo 

l> 

for every ^ S N for one (and hence for all, see [TH]) t G M. It is well-known that this 
initial value problem has unique global solutions which can be computed via the 

1^^^ ■ inverse scattering transform [16| . 

f~^ I The long-time asymptotics for this problem were first given by Novokshenov 

and HabibuUin [12 and were later made rigorous by Kamvissis [7], however, only 
in the case without solitons. The purpose of the present paper is to finally fill this 

^ ■ gap and show how to include solitons. As in [7], our approach is based on the 

JH I nonlinear steepest descent analysis for oscillatory Riemann-Hilbert problems from 

Deift and Zhou [1]. It turns out that in the case of solitons, two new phenomena 
enter the scene which require significant adaptions to the original method of Deift 
and Zhou. Of course our technique also applies to other soliton equations, e.g., the 
Korteweg-de Vries equation. 

First of all, it is well-known that there is a subtle nonuniqueness issue for the in- 
volved Riemann-Hilbert problems (see e.g. [H Chap. 38]). In fact, in certain excep- 
tional sets the corresponding vanishing Riemann-Hilbert problem has a nontrivial 
solution and hence by Fredholm theory the corresponding matrix Riemann-Hilbert 
problem has no solution at all. This problem does not affect the similarity region. 
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2 H. KRUGER AND G. TESCHL 

since it is easy to see that there it docs not happen for sufficiently large times. How- 
ever, in the soliton region, this occurs precisely in the neighborhoods of the single 
solitons. To avoid this problem we will work directly with the vector Riemann- 
Hilbert problem and impose a symmetry condition in order to ensure uniqueness. 
This also has the advantage that it eliminates the step of going forth and back 
between the vector and matrix Ricmann-Hilbert problem. It should be pointed 
out here that even the symmetry condition alone docs not guarantee uniqueness, 
rather we will need existence of a certain solution with an additional property. We 
will also demonstrate that this additional property is in fact necessary. 

Secondly, in the regions of the single solitons, the zeroth order asymptotics are 
not equal to zero but given by a one soliton solution. Hence for the usual perturba- 
tion argument based on the second resolvent identity to work, a uniform bound for 
the inverse of the singular integral equation associated with the one soliton solution 
is needed. Unfortunately, such a bound cannot be easily obtained. To overcome 
this problem we will shift the leading asymptotics from the one soliton solution 
to the inhomogcneous part of the singular integral equation and craft our Cauchy 
kernel in such a way that it preserves the pole conditions for this single soliton. 

To state our main result, we begin by recalling that the sequences a{n, t), h{n, t), 
n e Z, for fixed i G M, are uniquely determined by their scattering data, that 
is, by the right reflection coefficient i?+(z,i), \z\ = 1 and the eigenvalues Xj G 
(— cxD, —1) U (1, oo), j = 1, . . . , iV together with the corresponding right norming 
constants 7+j(f) > 0, j = 1,. . .,N. Rather than in the complex plane, we will 
work on the unit disc using the usual Joukowski transformation 

1 / 1 



(1.3) A=-(z + -j, z = \- VA2 - 1, AeC, |z|<l. 

In these new coordinates the eigenvalues Xj G (— cxd, — 1) U (1, oo) will be denoted by 
(j G (—1, 0) U (0, 1). The continuous spectrum [—1, 1] is mapped to the unit circle 
T. Moreover, the phase of the associated Riemann Hilbert problem is given by 

(1.4) $(z) = z-z"i + 2- log(z). 
and the stationary phase points, $'(z) = 0, are denoted by 



For J < — 1 we have zq G (0, 1), for — 1 < j < 1 we have zq G T (and hence 
Zq^ = "zo), and for j > 1 we have zo G (—1,0). For \j\ > 1 we will also need the 
value Co e (-1, 0) U (0, 1) defined via Re($(Co)) = 0, that is, 

(16) n^^Co-Q' 

^ ^ t 21og(|Co|) 

We will set Co = —1 if I jl < 1 for notational convenience. A simple analysis shows 
that for J < — 1 we have < Co < -^o < 1 and for j > 1 we have — 1 < zq < Co < 0. 
Furthermore, recall that the transmission coefficient T{z), \z\ < 1, is time inde- 
pendent and can be reconstructed using the Poisson- Jensen formula. In particular. 
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we define the partial transmission coefficient with respect to zq by 

r(z, zo) = 



(1.7) 



n lai^, 

ae(Co,o) 

n lai^ 

,ae(-i,o) 

n \C>^\'-7^]^-v(^J^ogi\Tis)\) 

,ae(-i,o)u(Co,i) / V T 



/ 

Zo 



expUl/log(in^)l)fSf h 



s-\-z ds 
s — z s 



20 e (-1,0), 
kol = 1, 
20 e (0,1). 



Here, in the case zq e T, the integral is to be taken along the arc ^{zq) — 
{z G T|Rc(z) < Rc(zo)} oriented counterclockwise. For zq g (—1,0) we set 
S(zo) = and for zq G (0,1) we set S(zo) = T. Then T{z,zq) is meromorphic 
for z £ C\I](zo). Observe T{z,zo) — T{z) once zq G (0,1) and (0,Co) contains no 
eigenvalues. Moreover, T(z, zq) can be computed in terms of the scattering data 
since \T{z)\^ = 1 - \R+{z,t)\^. 
Moreover, we set 



Toizo)^T{0,zo) 



(1.8) 



n lai 

ae(Co,o) 



20 e (-1,0), 



n iarMexpU./iog(|T(5)i)f 



|20 



and 



(1.9) 



^ae(-i,o) 

n lai 

k VCfce(-i,o)u(Co,i) 



Ti(zo) = — logT(z,zo) 
oz 

r E (c'-a) 

a6(Co,o) 

-1 



exp(^5i./log(|r(s)|)f^, zoe(0,l), 



j.J\og{\T{s)\y^ 



20 e (-1,0), 

|20| = 1, 



E^ ^ (C' - Ck) 

E iC^'^C,) + i-J\og{\T{s)\)^, zoe(0,l) 

^ae(-i,o)u(Co,i) T 



Theorem 1.1. Assume ()1.2p and abbreviate by Ck — — 2To~(T^T)' ^^^ velocity of 
the k'th soliton determined by Re($(Cfe)) = 0. Then the asymptotics in the soliton 
region, \n/t\ > 1 + C/ilog(i)^ for some C > 0, are as follows. 

Let e > sufficiently small such that the intervals [ck — £,Ck + s], 1 < A: < N, 
are disjoint and lie inside (— cxo, —1) U (1, cx)). 

If \j — Ck\ < £ for some k, one has 



Wi±7,(ML + 0(t-) 



(1.10) 






oo 



j=n+l 



2((7,(n,<)-l)C,^ + l) 



0(t-') 
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for any I > 1, where 

(1.11) 7fe(", t) = 7fcr(a, -cfc - ^/4~^r'^'^^'"^^'kl" 

If \j — Ck\ > £, for all k, one has 

CO 

l[{2a{j,t))^To{zo){l + 0{t-')), 

oo ^ 

(1.12) Y. b{j,t)^-T,{zo) + 0{t-'), 



j=n+l 



for any I > 1. 



In particular, we recover the well-known fact that the solution splits into a sum of 
independent solitons where the presence of the other solitons and the radiation part 
corresponding to the continuous spectrum manifests itself in phase shifts given by 
T(C/c, — Cfe — y/c^ — 1)~^. Indeed, notice that for (k 6 (—1,0) this term just contains 
product over the Blaschke factors corresponding to solitons Q with (k < Cj- For 
Ca; e (0, 1) we have the product over the Blaschke factors corresponding to solitons 
Q G (—1, 0), the integral over the full unit circle, plus the product over the Blaschke 
factors corresponding to solitons Q with C^ > Q . 

The proof will be given at the end of Section |4l Furthermore, in the remaining 
regions the analysis in Section 2] also shows that the Riemann-Hilbert problem re- 
duces to one without solitons. In fact, away from the soliton region, the asymptotics 
are given by 

oo oo 

l[{2a{j,t)) = To(-l) l[{2h{j,t)) (1 + 0{t-')) , 

oc ^ oc 

(1.13) Y. b{j,t)^^T,i-l)+ J2 b{j,t)+0{t-'), 

j=n+l j=n+l 

where d{n,t), b{n,t) are the solutions corresponding to the case without solitons 
and with _R+(z,0) replaced by 

(1.14) R+{z,0)^T{z,-l)-^R+{z,0). 
Note that the Blaschke product 

T(.,-i)= n la-i^ 

satisfies \T{z, — 1)| = 1 for z <E T. Hence everything is reduced to the case studied 
inH. 

Finally we remark that the same method can be used to handle solitons on a 
periodic background ^0\ (cf. also ^, [8], [9]). 

2. The Inverse scattering transform and the Riemann-Hilbert 

PROBLEM 

In this section we want to derive the Riemann-Hilbert problem from scattering 
theory. The special case without eigenvalues was first given in Kamvissis [7j. The 
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eigenvalues will be added by appropriate pole conditions which are then turned into 
jumps following Deift, Kamvissis, Kriecherbauer, and Zhou [5]. 

For the necessary results from scattering theory respectively the inverse scatter- 
ing transform for the Toda lattice we refer to [14] , [15] , [16] . 

Associated with a{t),b{t) is a self-adjoint Jacobi operator 

(2.1) H{t) = a{t)S+ + a- {t)S- + b{t) 

in i?^(Z), where S^f{n) — f^{n) = f{n± 1) are the usual shift operators and £^(Z) 
denotes the Hilbert space of square summable (complex-valued) sequences over Z. 
By our assumption (|1.2p the spectrum of H consists of an absolutely continuous part 
[—1, 1] plus a finite number of eigenvalues Afc G M\[— 1, 1], 1 < fc < A^. In addition, 
there exist two Jost functions 'ip±{z,n,t) which solve the recurrence equation 

(2.2) H{t)i^±{z,n,t)^ ^^^ i^±iz,n,t), \z\<l, 
and asymptotically look like the free solutions 

(2.3) lim z^"^P±{z,n,t)^l. 

n — *iboo 

Both ^±{z,n,t) are analytic for < |z| < 1 with smooth boundary values for 
\z\ = 1. The asymptotics of the two projections of the Jost function are 

(2.4) ^±iz,n,t) = ^^ (l + 2B±{n,t)z + 0{z^)), 
as z —^ 0, where 



(2.5) 



A+{n,t) = l[2a{j,t), B+{n,t) = - ^ b{j,t), 

n— 1 n— 1 

A-in,t)^ n MJ,t), B^{n,t) = - ^ 6(j,t). 

j=-oo j = -oo 

One has the scattering relations 



(2.6) T{z)ip^{z,n,t) ^ip±{z,n,t) + R±{z,t)^±{z,n,t), \z\ = I, 

where T{z), R±{z,t) are the transmission respectively reflection coefficients. The 
transmission and reflection coefficients have the following well-known properties: 

Lemma 2.1. The transmission coefficient T[z) has a meromorphic extension to 
the interior of the unit circle with simple poles at the images of the eigenvalues C,k ■ 
The residues ofT[z) are given by 

(2.7) Resc, T{z) = -Qk'^^^^ = -Ckl-At)Mt), 
where 

(2.8) 7±,fc(i)^'-El^±(C'=''^'*)l' 

nez 

and ip^{(k,n,t) ^ Hk{t)ip+iCk,n,t). 
Moreover, 



(2.9) T{z)R+{z,t) + T{z)R^{z,t) = 0, |T(z)|^ + |i?±(z, t)|^ = 1. 
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In particular one reflection coeflicient, say R{z,t) — R+{z,t), and one set of 
norming constants, say 7fc(t) — 7+,fe(t), sufflces. Moreover, the time dependence is 
given by: 

Lemma 2.2. The time evolutions of the quantities R+{z,t), 7+^fe(t) are given by 

(2.10) i?(z,i) = i?(z)e*(^-^"') 

(2.11) jk{t) ^ ^ke^^^^-^i^'^ 
where R{z) — R{z, 0) and 7^ ~ 7fe(0). 

Now we define the sectionally meromorphic vector 

^9 19^ .r,(.r,i\-f {T{z)^_{z,n,t)z^ ^+{z,n,t)z-^) , |z| < 1, 

{lAl) "^^^''^'^^-l (^,^(;2-i^„,i)2" r(z-i)^_(z-i,n,i)z-"), |z|>l. 

We are interested in the jump condition of m{z, n, t) on the unit circle T (ori- 
ented counterclockwise). To formulate our jump condition we use the following 
convention: When representing functions on T, the lower subscript denotes the 
non-tangential limit from different sides, 

(2.13) m±(z) = lim m(C), \z\ = 1. 

In general, for an oriented contour S, m^{z) (resp. m^{z)) will denote the limit 
of m{C) as C — > z from the positive (resp. negative) side of S. Using the notation 
above implicitly assumes that these limits exist in the sense that m{z) extends to 
a continuous function on the boundary. 

Theorem 2.3 (Vector Riemann-Hilbert problem). Let S+{H{{))) = {R{z), \z\ — 

1; (Cfe:7fc), 1 < fc < N} the right scattering data of the operator H{0). Then 
m(z) — m(z,n,t) defined in (|2.12[) is meromorphic away from the unit circle with 
simple poles at Ck, Ck '^"■'^ satisfies: 
(i) The jump condition 



(2.14) m+yz) — m-[z)vyz), v{z) = ' 1 \ 'i y ' 



i?(z)e**(^) 1 



for z € T, 
(ii) the pole conditions 



(2.15) 



-<""'='^.^."'=' -a.:-. 



Res,- 

'■fc 



m,{z) — lim m,{z) '^ 

iii) the symmetry condition 



(2.16) m(z-i) = m(z) (^J J^ 

(iv) and the normalization 

(2.17) 771(0) = {mi 7712), Til • ^2 = 1 rni > 0. 
Here the phase is given by 

(2.18) $(z) = z-z"^ + 2-log z. 
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Proof. The jump condition (|2.14p is a simple calculation using the scattering re- 
lations (|2.6|) plus ()2.9|1 . The pole conditions follow since T{z) is meromorphic in 
I z I < 1 with simple poles at (k and residues given by (|2.7p . The symmetry condition 
holds by construction and the normalization (|2.17p is immediate from the following 
lemma. D 

Observe that the pole condition at Cfc is sufficient since the one at (^^ follows 
by symmetry. 

Moreover, we have the following asymptotic behaviour near z = 0: 

Lemma 2.4. The function m{z,n,t) defined in p.l2|) satisfies 

(2.19) miz,n,t)=(^A{n,t){l-2B{n-l,t)z) -^(1 + 2B(n,i)z)) + 0(^2). 

Here A{n,t) = A^{n,t) and B(n,t) — B^{n,t) are defined in (12. 5p . 

Proof ThisfoUowsfromdlH) andr(z) = A+A_(l-2(B+-fe + B_)z-|-0(z2)). Q 

For our further analysis it will be convenient to rewrite the pole condition as a 
jump condition and hence turn our meromorphic Riemann-Hilbert problem into a 
holomorphic Riemann-Hilbert problem following [5]. Choose e so small that the 
discs |z — Cfcl < e are inside the unit circle and do not intersect. Then redefine m 
in a neighborhood of (k respectively (^^ according to 



™W f.^.e'^^-c) , ' \z-Ck\<e, 



(2.20) m{z) = < 




m(z) -C.7 1 |z-i-Cfc|<e, 



m{z), else. 

Then a straightforward calculation using Res,;; m = \imz—,({z — ()m{z) shows 

Lemma 2.5. Suppose m{z) is redefined as in (|2.20|) . Then m{z) is holomorphic 
away from the unit circle and satisfies (I2.14p . ()2.16|) . ()2.17|) and the pole conditions 
are replaced by the jump conditions 

0\ 




(2.21) 

TO+(z) = m_(z) I " ^-C,:' I, \z-^-Ck\=£, 

where the small circle around C,k is oriented counterclockwise and the one around 
C^ is oriented clockwise. 

Next we turn to uniqueness of the solution of this vector Riemann-Hilbert prob- 
lem. This will also explain the reason for our symmetry condition. We begin 
by observing that if there is a point zi G C, such that m{zi) = (O O), then 
^(•2^) = z-z '^{z) satisfies the same jump and pole conditions as m,{z). However, 
it will clearly violate the symmetry condition! Hence, without the symmetry con- 
dition, the solution of our vector Riemann-Hilbert problem will not be unique in 
such a situation. Moreover, a look at the one soliton solution verifies that this case 
indeed can happen. 
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Lemma 2.6 (One soliton solution). Suppose there is only one eigenvalue and a 
vanishing reflection coefficient, that is, S+{H{t)) = {R{z) = 0, \z\ — 1; (Ct)} with 
C e (-1,0) U (0,1) and 7 > 0. Then the Riemann-Hilbert problem ^J^-^JJii 
has a unique solution is given by 

(2.22) 

mo(z)-(/(z) fil/z)) 

f{z) = _\ (^7(71, i)C'^^^ + 1 - C' 

^i-e+^[n,t)^i-e+ei{n,t) V'' ' ^-c 

where ^{n,t) = 76*^^'. /n particular, 



..„^ 4 r .^ l-C2+7(n,^) 7(n,t)C(C'-l) 

(2.23) A+(n,i) = W- — - — - — — -, B+{n,t)- 



l-C2+7(n,t)C2' ^^ ' ' 2{\-e^l{n,t)ey 

Proof. By symmetry, the solution must be of the form mQ{z) = (f{z) f{\/z)\ 
where f{z) is meromorphic in C U {00} with the only possible pole at C,. Hence 



A\ z-C 

where the unknown constants A and B are uniquely determined by the pole con- 
dition Res^ /(z) = —C,"l{n,t)f{(^~'^) and the normalization /(0)/(oo) — 1, /(O) > 
0. D 

In fact, observe f{zi) = f{zf^) = if and only if zi = ±1 and 7 = ±(C "'^ — C)- 
Furthermore, even in the general case m{zi) = (O O) can only occur at zi = ±1 
as the following lemma shows. 

Lemma 2.7. //m(zi) — (O O) for m defined as in (|2.12p . then zi = ±1. More- 
over, the zero of at least one component is simple in this case. 

Proof. By (|2.12p the condition m{zi) — (O O) implies that either the Jost solutions 
tp-izi, n) and ?/'+(zi, n) are linearly dependent or T{zi) = 0. This can only happen, 
at a band edge, zi = ±1, or at an eigenvalue zi = Q. 

We begin with the case zi = Q. In this case the derivative of the Wron- 
skian W{z) = a{n){ip+{z,n)^-{z,n + 1) — ^+[z,n + l)^_(z,n)) does not van- 
ish ^W^(2^)|z=zi 7^ (Uni Chap. 10]). Moreover, the diagonal Green's function 
g{X, n) = W{z)^^'ip^{z, n)xj}_{z, n) is Herglotz and hence can have at most a simple 
zero at z = zi. Hence, if '(/'+(Cj, '^) — ^-(Cj, '^) — 0, both can have at most a simple 
zero at z = Q. But T{z) has a simple pole at Q and hence T {z)i}) -{z , n) cannot 
vanish at z = ^j, a contradiction. 

It remains to show that one zero is simple in the case zi = ±1. In fact, one 
can show that j^W{z)\z=zi ^ in this case as follows: First of all note that V'±(^) 
(where / denotes the derivative with respect to z) again solves Hijj'^{zi) = AiV'±(zi) 
if zi = ±1. Moreover, by W{zi) = we have V'-i-(^i) = ciIj-{zi) for some constant 
c (independent of n) . Thus we can compute 

I^'(zi) = I^(^^(zi),^_(zi)) + T4^(^+(zi),^L(zi)) 

= c-iM/(V;;(zi),^+(zi))+cTy(^_(zi),^^(zi)) 
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by letting n — > +00 for the first and n — > — cxd for the second Wronskian (in which 
case we can replace '4'±{zi) by z-^ "), which gives 

W'iz.) = ^. 

Hence the Wronskian has a simple zero. But if both functions had more than simple 
zeros, so would the Wronskian, a contradiction. D 

Finally, it is interesting to note that the assumptions C, G ( — 1,0) U (0,1) and 
7 > are crucial for uniqueness. Indeed, if we choose 7 = C^ — 1 < 0, then every 
solution is a multiple of f{z) = zC,~^(z — C,)~^ which cannot be normalized at 0. 

3. A UNIQUENESS RESULT FOR SYMMETRIC VECTOR RiEMANN-HlLBERT 

PROBLEMS 

In this section we want to investigate uniqueness for the holomorphic vector 
Riemann-Hilbert problem 

mj^{z) — m_(z)u(z), z G S, 

(3.1) m(z"i) = m(z)L ^ 

m(0) = (1 TO2) . 

where S is a nice oriented contour (see Hypothesis I A. 1 p . symmetric with respect to 
z 1-^ z~^ , and v is continuous satisfying det(w) ~ 1 and 

(3.2) v{z-^)={^^ J)^(^)"'(! J)' ^^^- 

The normalization used here will be more convenient than (|2.17p . In fact, (|2.17p 
will be satisfied by ttt-j m{z). 

Now we are ready to show that the symmetry condition in fact guarantees unique- 
ness. 

Theorem 3.1. Suppose there exists a solution m{z) of the Riemann-Hilbert prob- 
lem p.ip for which m{z) — (0 O) can happen at most for z = ±1 in which case 
limsup2_,_l_]^ TwIT *'^ bounded from any direction for j — 1 or j — 2. 

Then the Riemann-Hilbert problem (|3.ip with norming condition replaced by 

(3.3) to(0) = (a 7712) 

for given a G C, has a unique solution ma{z) = am(z). 

Proof. Let ma{z) be a solution of p.ip normalized according to (|3.3p . Then we 
can construct a matrix valued solution via M = [m, ma) and there are two possible 
cases: Either det M{z) is nonzero for some z or it vanishes identically. 

We start with the first case. Since the determinant of our Riemann-Hilbert prob- 
lem has no jump and is bounded at infinity, it is constant. But taking determinants 
in 



M(z-i) = Af (z) (^5 J) 



gives a contradiction. 

It remains to investigate the case where det(Af) = 0. In this case we have 
Tna{z) = d{z)m{z) with a scalar function S. Moreover, S{z) must be holomorphic 
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for z e C\I] and continuous for z G E except possibly at the points where to(zo) = 
(O O). Since it has no jump across E, 

6+{z)m^{z) — mQ_+(z) — ■ma,^{z)v{z) — S-{z)m^{z)v{z) = (5_(z)m+(z), 

it is even holomorphic in C\{±1} with at most simple poles at z = ±1. Hence it 
must be of the form 

c-/ \ . B C 

z — 1 z + 1 
Since <5 has to be symmetric, S{z) = 6{z~^), we obtain B = C = 0. Now, by the 
normalization we obtain d{z) — A = a. This finishes the proof. D 

Furthermore, note that the requirements cannot be relaxed to allow (e.g.) second 
order zeros in stead of simple zeros. In fact, if r7T,(z) is a solution for which both 
components vanish of second order at, say, z = +1, then m{z) = , _f,-,2 ^7i(z) is a 
nontrivial symmetric solution of the vanishing problem (i.e. for a = 0). 

By Lemma 12.71 we have 

Corollary 3.2. The function m{z, n, t) defined in (j2.12p is the only solution of the 
vector Riemann-Hilbert problem (J2.14p - (j2.17p . 

Observe that there is nothing special about the point z = where we normalize, 
any other point would do as well. However, observe that for the one soliton solution 
(HSU, f{z) vanishes at 

.- 1 + 7-C^ 

' ^ (7 - i)e + 1 

and hence the Riemann-Hilbert problem normalized at this point has a nontrivial 
solution for a — and hence, by our uniqueness result, no solution for a = I. This 
shows that uniqueness and existence connected, a fact which is not surprising since 
our Riemann-Hilbert problem is equivalent to a singular integral equation which is 
Fredholm of index zero (see Appendix E| • 

4. SOLITONS AND THE SOLITON REGION 

This section demonstrates the basic method of passing from a Riemann-Hilbert 
problem involving solitons to one without. Furthermore, the asymptotics inside 
the soliton region are computed. Solitons are represented in a Riemann-Hilbert 
problem by pole conditions, for this reason we will further study how poles can be 
dealt with in this section. 

For easy reference we note the following result which can be checked by a straight- 
forward calculation. 

Lemma 4.1 (Conjugation). Assume that E C E. Let D be a matrix of the form 

U{z)-^ 



(4-1) D[z)-^ ^ ^(^, 

where d : C\E —> C- is a sectionally analytic function. Set 

(4.2) to(z) = vi{z)D{z), 
then the jump matrix transforms according to 

(4.3) v{z) = D^{zy^v{z)D+{z). 
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// d satisfies d{z^^) = d(z)^^ and d{0) > 0. Then the transformation rh{z) = 
m{z)D(z) respects our symmetry, that is, rh{z) satisfies (|2.16|) if and only if m{z) 
does. 



In particular, we obtain 

(4.4) u = 
respectively 

(4.5) V 



Wll Vl2d^ 

V2id~^ '(;22 



d^vii vi2d+d- 

V2ld+ d_ d^«22 



z e E\S, 



ze snE. 



In order to remove the poles there are two cases to distinguish. If Afc > ^ (Co+CcT^) 
the jump is exponentially decaying and there is nothing to do. 

Otherwise we use conjugation to turn the jumps into exponentially decaying 
ones, again following Deift, Kamvissis, Kriecherbauer, and Zhou 5J. It turns out 
that we will have to handle the poles at C^ and C,^ in one step in order to preserve 
symmetry and in order to not add additional poles elsewhere. 

Lemma 4.2. Assume that the Riemann-Hilbert problem for m has jump conditions 
near ^ and i^^^ given by 



(4.6) 



1 



m+{z) ^m-{z) _^ 
\z~c 

m^{z) = ?7i_(z) 







1 



k-Cl=£, 



-'-C\=e. 



Then this Riemann-Hilbert problem is equivalent to a Riemann-Hilbert problem for 
fa which has jump conditions near C, and C,^^ given by 



Th+[z) = rh-{z) 







C(z-C)7 



1 



m+{z)=m^{z)\ (^_g)2 , 

\Cz(Cz-l)l ^1 

and all remaining data conjugated (as in Lemma \4.1\ l by 
(4.7) D{z) ^ ( ^ ^ j . 



\z-'-C\^e, 







Proof. To turn 7 into 7 ^, introduce D by 



Diz) 




and note that D{z) is analytic away from the two circles. Now set rh{z) = 
m{z)D{z), which is again symmetric by D{z^^) — {io)D{z){iq). The jumps 
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along |z — CI — ^ ^-nd \z^^ — C\ = e follow by a straightforward calculatfon and the 
remaining jumps follow from Lemma [4711 D 

Now we are ready to prove our main result: 

Proof of Theorem 11. Jl We begin by observing that the partial transmission co- 
efficient T{z, zq) introduced in (|1.7p satisfies the following scalar meromorphic 
Riemann-Hilbert problem: 

(i) T{z,zq) is meromorphic in C\S(zo), where S(zo) is the arc given by 
T,{zq) = {z e T|Re(z) < Re(zo)}, with simple poles at Cfc and simple 
zeros at Q^^ for all k with Afc < ^(Co + CcT^); 
(ii) r+(z,zo) = r_(z, zo)(l - |i?(z)|2) for z € S(zo), 
(iii) r(z-i,zo) = r(z,zo)-\ z e C\I](zo), and T(0, zq) > 0. 
Note also r(z, zq) — r(z, zq) and in particular T{z, zq) is real- valued for z G M. 
Next introduce 

z-Cfc 



D{z) 



where 



1 

_ Cfc7fce_ 



-(c.) ^'"^'=1 " ^o(^), |z - al < e, Afc < i(Co + Co"'), 



z-Cfc 






, , ^oW, |z-i-Cfc|<£, Afe<i(Co+Co"'), 

Do(z), else, 



Note that we have 






Now we conjugate our vector m(z) defined in (|2.12p respectively (I2.20p using D{z), 
rh{z) = m{z)D{z). Since T{z, zq) is either nonzero and continuous near z = ±1 (if 
±1 ^ S(zo)) or it has the same behaviour as T{z) near z = ±1 (if ±1 G S(zo)), the 
new vector rh{z) is again continuous near z = ±1 (even if T{z) vanishes there). 

Then using Lemma 1411 and Lemma 1421 the jump corresponding Afc < ^(Co + Ccr') 
(if any) is given by 



(4.8) 



v{z) = ^ a7.T(.,.o^^e-(c.) 1 \z-a\^e, 



and corresponding Afc > i(Co + Co ') (i^ any) by 

1 0\ 

a-7.T(.,^o)~^c'-c:>=) J, |z-Cfc|=e, 

(4 9) ' ^~''' ^ 

In particular, an investigation of the sign of Re(<i>(z)) shows that all off-diagonal 
entries of these jump matrices, except for possibly one if Cfeo — Co ior some fco) are 
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exponentially decreasing. In this case we will keep the pole condition which now 
reads 

Resr^ m(z) = lim m(z) ( , ,_2 tis>(Ci ) n ) ' 



(4.10) 



Res,-i to(z) = lim rh(z) „ " 

'■'=0 ^^(-^ ' ' \0 



Ck'ikoT{Cko,zor\"'^'^^'>) 



Furthermore, the jump along T is given by 

..... -.. , /6_(z)-i6+(z), A(z)>Ao, 

where 






and 




Now the jump along T can also be made arbitrarily small following the nonlinear 
steepest descent method developed by Deift and Zhou Q: Split the Fourier trans- 
form of the reflection coefficient into a part which has an analytic extension to a 
neighborhood of the unit circle plus a small error. Move the analytic part away 
from the unit circle using the factorizations from above. Since the Fourier coeffi- 
cients decay faster than any polynomial, the errors from both parts can be made 
0(i~') for any I G N. We refer to [4j respectively [7] for details. Hence we can 
apply Theorem lA. 61 as follows: 

If I J — Cfcl > e for all k we can choose 70 = and Wq =0. Since the error between 
w* and Wq is exponentially small, this proves the second part of Theorem II .11 upon 
comparing 

(4.12) ^(.).^(.)(^(^^^°) ^(^;^^^_ 

with dmn). 

Otherwise, if | j — Cfe| < e for some fc, we choose 70 = 7fe(f^, t) and Wq = 0. Again 
we conclude that the error between w* and Wq is exponentially small, proving the 
first part of Theorem 11.11 D 
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Appendix A. Singular integral equations 

In this section we show how to transform a meromorphic vector Riemann-Hilbert 
problem with simple poles at (, (~^, 



m+{z) = m_(z)u(z), 


z e E, 


(A.l) 




Res^ to(z) — lim m(z) 


0' 

"C7 0, 


m(z"^)=m(z)L A 


; 


m(0) = (1 ma) , 





ReS(;-i m{z) = lim m{z) 



rw 





where ^ € (—1,0) U (0, 1) and 7 > 0, into a singular integral equation. Since we 
require the symmetry condition (|2.16p for our Riemann-Hilbert problems we need 
to adapt the usual Cauchy kernel to preserve this symmetry. Moreover, we keep 
the single soliton as an inhomogeneous term which will play the role of the leading 
asymptotics in our applications. 

Hypothesis H. A.l. Let E consist of a finite number of smooth oriented finite 
curves in C which intersect at most finitely many times with all intersections being 
transversal. Assume that the contour E does not contain 0, (^ and is invariant 
under z >— > z'~^ . It is oriented such that under the mapping z i— > z^^ sequences 
converging from the positive sided to S are mapped to sequences converging to the 
negative side. Moreover, suppose the jump matrix v can be factorized according to 
V — bZ b+ = {I — W-)~^{I + Wji-), where w± — ±(6± —I) are continuous and satisfy 

(A.2) w±{z-')=(^^ J)^t(^)(? nV ^eS. 







The classical Cauchy-transform of a function / : S ^ C which is square inte- 
grable is the analytic function Cf : C\S] -^ C given by 

(A.3) Cf{z) ^^ f I^ds, z e C\S. 

Denote the non-tangential boundary values from both sides (taken possibly in the 
I/^-sense — see e.g. [21 eq. (7.2)]) by C+/ respectively C_/. Then it is well-known 
that C+ and C_ are bounded operators i^(S) -^ L^(E), which satisfy C+ — C_ — I 
and Cj^C- = (see e.g. [T]). Moreover, one has the Plemelj-Sokhotsky formula 

C± = \{\H±l), 
where 

(A.4) Hf{t) - -j l^ds, t G E, 

ttJy t - s 

is the Hilbert transform and j^ denotes the principal value integral. 

In order to respect the symmetry condition we will restrict our attention to the 
set i^(E) of square integrable functions / : E -^ C^ such that 

(A.5) /(z-i) = /(z) (J I 
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Clearly this will only be possible if we require our jump data to be symmetric as 
well (i.e., Hypothesis I A. II holds) . 

Next we introduce the Cauchy operator 

(A.6) (Cf){z)^^ f f{s)n^{s,z) 

acting on vector- valued functions / ; S ^ C^. Here the Cauchy kernel is given by 

/ z-c' 1 n \ /^ ^— \ 

(A.7) 17c(.,z)= h-C--^ -C_^r^=r"' 0'~^" ^ ^r^' 

for some fixed C, ^Ti. In the case C = cxo we set 

1 1 



(A.8) r!oo(s,z)= ( ^--Q ^ ^jds. 

and one easily checks the symmetry property: 

(A.9) fic(lMlA)=(j J)«c(^:^)(j J)- 

The properties of C are summarized in the next lemma. 

Lemma A. 2. Assume Hvpothesis \A.l[ The Cauchy operator C has the properties, 
that the boundary values C± are bounded operators -^^(E) ^> ^^(S) which satisfy 

(A. 10) C+ - C- = I 

and 

(A.ll) (C/)(ri) = (0 *), (C/)(C)-(* 0). 

Furthermore, C restricts to L^(E), that is 

(A.12) (C/)(z-i) = (a/)(z)(^J J), ^eC\E 

/or / e is(S) and if w± satisfy fH \A.l\) we also have 

(A.13) C±{fw^){l/z)^C^{fw±)iz)(^^ J^, zeE. 

Proof. Everything follows from (jA.Op and the fact that C inherits all properties 
from the classical Cauchy operator. D 

We have thus obtained a Cauchy transform with the required properties. Fol- 
lowing Section 7 and 8 of [T] , we can solve our Riemann-Hilbert problem using this 
Cauchy operator. 

Introduce the operator C„, : ^^(E) -^ L'^i^) by 

(A.14) C^f^C+{fw-) + C-ifw+), /eL2(S) 

and recall from Lemma l2.6l that the unique solution corresponding to w = I is given 

by 

moiz) = {f{z) /(i)), fiz) = ^-^(^^C'"-^ + l-e) 

Observe that for 7 = we have f{z) = 1 and for 7 = 00 we have f{z) = C,'^ ^^f . . 
In particular, too(z) is uniformly bounded away from C, for all 7 G [0, 00]. 
Then we have the next result. 
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Theorem A. 3. Assume Hypothesis \ A. 1[ 

Suppose m solves the Riemann-Hilbert problem ()A.1|1 . Then 

(A. 15) m(2;) = (l-co)mo(z) + — ^ / ^i{s){w+{s) + W-{s))VLq{s,z), 

where 

/i = TO+&7 = m^bZ and Cq — [ t;-t / A'('S)(w+(s) + i(;_(s))rir(s, 0) I . 

V27riJs /i 

Here {m)j denotes the j 'th component of a vector. Furthermore, /i solves 
(A.16) (I-C^)Ai=(l-co)mo(z). 

Conversely, suppose jl solves 
(A.17) (I - C^)/i = mo(z), 

anrf 

5o- (^^^M(s)(u;+(s)+ii;_(,s))r!cGs,0)) ^-1, 

then m defined via (jA.lSp . wit/i (1 — cq) = (1 + cq)^^ and /i = (1 + cq)^"'^/!, solves 
the Riemann-Hilbert problem (jA.ip arirf /i = m±b^ . 

Proof. This can be shown as in the non-symmetric case (cf. e.g. [3]). D 

Note that in the special case 7 = we have mQ{z) = (l l) and we can choose 
C as we please, say ^ = 00 such that cq ~ cq — in the above theorem. 

Hence we have a formula for the solution of our Riemann-Hilbert problem m{z) 
in terms of {I—Cw)~^'mo and this clearly raises the question of bounded invertibility 
of I — Cw This follows from Fredholm theory (cf. e.g. [19 J: 

Lemma A. 4. Assume Hypothesis \ A. 1\ The operator I — Cw is Fredholm of index 
zero, 

(A. 18) ind(I - Cu,) = 0. 

By the Fredholm alternative, it follows that to show the bounded invertibility of 
I — Cw we only need to show that ker(I — Cw) — 0. The latter being equivalent to 
unique solvability of the corresponding vanishing Riemann-Hilbert problem in the 
case 7 = (where we can choose ^ = 00 such that Cq = Cq = 0). 

Corollary A. 5. Assume Hypothesis \A.1[ A unique solution of the Riemann- 
Hilbert problem (jA.ip with 7 = exists if and only if the corresponding vanish- 
ing Riemann-Hilbert problem,, where the normalization condition is replaced by 
m(0) = (0 TO2), with TO2 arbitrary, has at most one solution. 

We are interested in comparing two Riemann-Hilbert problems associated with 
respective jumps wo and w with \\w — wo|!oo small, where 

(A.19) ||u;||co - Ik+||L~(S) + Ilii'-llL-CS)- 

For such a situation we have the following result: 

Theorem A. 6. Assume that for some data (q, Jq, vJq the operator 

(A.20) I-C^* :L2(s)->l2(S) 

has a bounded inverse, where the bound is independent of t, and suppose that the 
corresponding Cg q is away from —1 again uniformly in t. 
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Furthermore, let C = Co; 7* = 7o ^'^'^ assume w* satisfies 

(A.21) |k*-wolloo <a(t) 

for some function a{t) — > as i — » oo. Then (I — C^')"^ ■ L'j.{Yi) -^ -^s(^) '^^^'^ 
exists for sufficiently large t and the associated solutions of the Riemann-Hilbert 
problems (jA.ip only differ by 0{a{t)) uniformly in z away from E U {C, C""'^}- 

Proof. By boundedness of the Cauchy transform, one has 

\\{C^,t - C^tJW < C07isi||w||oo- 

Thus, by the second resolvent identity, we infer that (I — C^*)'^'^ exists for large t 
and 

||(i-c^O"'-(i-c^^)"'ll-o(a(t)). 

From which the claim follows since this implies \cq — Cq q\ — 0{a{t)) respectively 
IIm* — Molli^ = 0{a{t)) and thus TO*(2:) — mo(z) = 0{a{t)) uniformly in away from 

su{c,r'}- □ 
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